Ferromagnetic transition in a double-exchange system containing impurities 



O 

o 

(N 



CO 

c3 



I 

o 
o 



(N 
> 
O 

in 
in 
in 
o 



c3 



O 

o 



X 
S3 



Mark Auslender 1 and Eugene Kogan 2 

1 Department of Electrical and Computer Engineering, Ben-Gurion University of the Negev, P.O.B. 653, Beer-Sheva, 84105 

Israel 

2 Jack and Pearl Resnick Institute of Advanced Techno 



Department of Physics, Bar-Ilan University, Ramat-Gan 52900, 
Israel 
(February 1, 2008) 

We study ferromagnetic transition in three-dimensional double-exchange model containing im- 
purities. The influence of both spin fluctuations and impurity potential on conduction electrons is 
described in coherent potential approximation. In the framework of thermodynamic approach we 
construct Landau functional for the system "electrons (in disordered environment) + core spins". 
Analyzing the Landau functional we calculate the temperature of ferromagnetic transition Tc and 
paramagnetic susceptibility x- F° r To, we thus extend the result obtained by Furukawa in the 
framework of the Dynamical Mean Field Approximation, with which our result coincides in the 
limit of zero impurity potential. We find, that the alloy disorder, able to produce a gap in density 
of electron states, can substantially decrease Tc with respect to the case of no impurities. 

PACS numbers: 75.10.Hk, 75.30.Mb, 75.30. Vn 



I. INTRODUCTION 

The recent rediscovery of colossal magnetoresis- 
tance in doped Mn oxides with perovskite structure 
Ri-zDzMnOa (R is a rare-earth metal and D is a di- 
valent metal, typically Ba, Sr or Ca) substantially in- 
creased interest in the double-exchange (DE) model ||[| . 
Several approaches were used lately to study the thermo- 
dynamic properties of the DE model, including the Dy- 
namical Mean Field Approximation (DMFA) j| (and ref- 
erences therein; DMFA itself see [pi), Green functions de- 
coupling techniques 0, Schwinger bosons Q| , variational 
mean- field approach ]8| and numerical methods [p|-pd|. 
In all these approaches chemical disorder introduced by 
doping, which is generic for the manganites, has not been 
taken into account. 

Recently we have shown that the concurrent action of 
the chemical and magnetic disorder, is crucial for the the 
description of the density of states and conductivity in 
manganites |l2| . In the present paper we consider the fer- 
romagnetic transition in the case of non-zero potential of 
randomly distributed impurities, using the same coherent 



as m our previ- 
Briefly, the 



potential approximation (CPA) |13| 
ous paper Q (see also relevant Ref. 
effect of impurities on the transition is the reduction of 
Tc as the impurity potential strength increases. In more 
detail, at values of the impurity potential at which the 
electron band is split off into two sub-bands, the mecha- 
nism of ferromagnetic exchange is other than in the case 
of the zero or weak potential. 



II. HAMILTONIAN AND CPA EQUATIONS 

We consider the DE model with the inclusion of the 
single-site impurity potential. We apply the quasiclassi- 



cal adiabatic approximation and consider each core spin 
as a static vector of fixed length S (Sj = where rij 
is a unit vector). The Hamiltonian of the model if([rij]) 
in site representation is 



Hi 



U-j + (ej - Jn, ■ a) 5i 



(1) 



where ti-j is the electron hopping, q is the on-site en- 
ergy, J is the effective exchange coupling between a core 
spin and a conduction electron and a is the vector of the 
Pauli matrices. The hat above the operator reminds that 
it is a 2 x 2 matrix in the spin space (we discard the hat 
when the operator is a scalar matrix in the spin space). 
The Hamiltonian (^) is random due to randomness of a 
core spin configuration [rij] and the randomness of the 
on-site energies e,. 

To handle CPA we present Hamiltonian (Q) as 



H" 



H 



(2) 



(the site independent self-energy S(-E') is to be deter- 
mined later), and construct a perturbation theory with 
respect to random potential V. To do this let us intro- 
duce the exact T-matrix as the solution of the equation 



T = V + VGnT, 



(3) 



in which Go = (E — H°) 1 . For the exact Green function 
G we get 

G = G + G TG . (4) 
The self energy is determined from the requirement 

'g) = Go (5) 



1 



where the angular brackets with the indexes mean av- 
eraging over the configurations of both core spins and 
impurities. In CPA T is considered in a single-site ap- 
proximation, at which Eq. (^) is reduced to the equation 



Ti) = 0, 



(6) 

where Tj is the solution of the equation 

% = Vi + V t g{E)f t . (7) 

Here 

g(E)=g (E-t), (8) 

(9) 



g (E) = I — de, 



E-s 



Nq{e) being the bare (i.e. for e, = and J = 0) density 
of states (DOS). Finally Eq. (||) can be transformed to 
an algebraic equation for the 2x2 matrix E: 



g = (lV + §- l -Vi 



III. BAND STRUCTURE 

To consider the evolution of the DOS 
1 . 



(10) 



N(E) = -Im 5c , 



(11) 



where g c = Tr g (here Tr means the trace over spin states 
only), with the variation of the impurity concentration 
and potential strength, we exploit the semi-circular (SC) 
bare DOS given at |e| < W (W is half of the bandwidth) 
by 



Ar ( £ ) = J_Ji-(±) 



(12) 



and equal to zero otherwise. For this DOS (which is exact 
on a Caley tree) 



9o(E) = - 




(13) 



Hence we obtain 

t = E-2wg-g-\ (14) 
where w = W 2 /8. Thus, Eq. ( |l0| ) transforms to 

g= ((E-e+ Jna-2wg)~ 1 ) . (15) 

\ / n,e 

It is convenient to write the locator g in the form 



(16) 



where I is a unity matrix. For the charge locator g c and 
spin locator g s we obtain the system of equations 



E — e — wg c 



(E — e — wg c ) — ( Jn — w. 



wg s - Jn 



(E-e- wg c ) - (Jn - wg s ) 2 




(17) 



In the strong Hund coupling limit (J — » oo) we obtain 
from Eqs. (|l^) two decoupled spin sub-bands. For the 
lower sub-band, after shifting the energy by J we obtain 



1 



E — e — wg c — wn z g 
n z 

E — e — wg c — wn z g s 



(18) 



where g s = (0,0, g s ) an axis OZ is directed along the 
average magnetization of core spins m) . 

In fact, details of alloying define how to average over 
the configurations of impurities in Eqs. ( p7| , ^8| ). We use 
for this random substitution model of disorder. That 
is ti — V with the probability x and e$ = with the 
probability 1 — x, where V and x are the impurity po- 
tential and concentration, respectively. As to core spins, 
once CPA is introduced, their configuration probability 
should be determined self-consistently in order to close 
Eqs. ( |l7| , p~8| ) . We have proved elsewhere that in the pres- 
ence of an annealed disorder (including dynamical one) 
DMFA ansatz for the disorder configuration probability 
|| keeps the free energy stationary against variations of 
CPAE@. 

For the following two cases, however, only the averag- 
ing over e is left. For a saturated ferromagnetic (FM) 
phase (m = 1), we obtain g s — g Cl and closed equation 
for the charge locator is 



9c 



2wg c 



(19) 



For a paramagnetic (PM) phase (m = 0), we obtain 
g s = 0, and closed equation for the charge locator is 



7 (o) _ 



1 



E — e — wg c 



(0) 



(20) 



It appears that calculation of gi°^ is sufficient for obtain- 
ing Tq and the PM susceptibility x as functions of x and 
V. 

The results of numerical calculation of the DOS at the 
PM state are presented on Fig. 1 (specific x and different 
V). 
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FIG. 1. The calculated DOS as function of the relative 
strength of impurity potential V/W for x = 0.18. 

and on Fig. 2 (specific V and different x). 



The dependence of the gap A on x at several V/W is 
shown on on Fig. 3, which summarizes the regularities 
discussed above. 
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FIG. 3. The calculated gap in the DOS as function of im- 
purity concentration x at different strengths of impurity po- 
tential V. 




FIG. 2. The calculated DOS as function of impurity con- 
centration x for V/W = 0.69. 



When comparing Eqs. (|19|) and (|20|) it is seen that 
the DOS in the FM phase is equal to the DOS in the 
PM state of another model, with increased by a factor 
of y/2 bare bandwidth. Using this property and Fig. 1 
we conclude that at an appropriate V/W there may be a 
gap between the conduction band states and the impu- 
rity band states in the PM phase, while the FM DOS is 
gapless. This may explain metal-insulator transition ob- 
served in manganites and magnetic semiconductors [p"2[ . 
It is also seen from Fig. 2 that the gap in the DOS exist- 
ing at low concentrations may close at higher concentra- 
tions. The value of the impurity potential V c = W/\[2 
detaches two types of the gap behavior. At V < V c the 
gap opens for some x < 0.5 or does not open at all, at 
V > V c the gap exists for all < x < 0.5, and at V = V c 
the gap closes exactly for x = 0.5. Anyhow, even if the 
gap is closed there still may exist pseudogap at a strong 
enough potential (see Fig. 2). 



It is worth noting, that at V > V c , A is near to saturate 
after x > 0.25. This may explain the near independence 
of the resistivity activation energy on x in the PM phase 
observed in single crystalline manganites (the fact 
used by some authors to support the polaron scenarios 
of the transport at T > Tc). 



IV. LANDAU FUNCTIONAL 

In our approach the calculation of thermodynamic 
properties is based on the analysis of the Landau func- 
tional. We start from the exact partition function of the 
system electrons + core spins 



exp 



-/3£l([n i ])+/JH.£;i* 



j=i 



;!</» (2i) 



where N is the number of core spins, H is magnetic field, 
j3 = 1/T, and f2([iij]) is the grand canonical potential of 
the electron subsystem for a given core spin configura- 
tion [rij]. Our aim is to obtain Landau functional of the 
system. 

Since we use CPA for electronic properties it is consis- 
tent to construct the Landau functional using a mean- 
field approach. Mean-field means that we do not take 
into account large-scale fluctuations of the macroscopic 
magnetization M = Y2j n j > but within CPA we cer- 
tainly take into account microscopic fluctuations of rv 
In the mean-field approximation all energy levels depend 
only upon M. Hence we may approximately put 



fi([ nj ])«fi(Af), 



(22) 



3 



where f2(m) is the grand canonical potential calculated 
within CPA at a non-zero m. Then the partition function 
can be written as 

exp [-/3L(M, H)] dM, (23) 
where 

L(M, H) = n(M) - iVH • M - T5(M); (24) 
the quantity 

AT \ AT 
i=l / »=1 



S(M) = In / <5 



is the entropy of the core spin subsystem. So we may 
identify the functional L(M, H) in the exponent of Eq . 
( f23| ) with the Landau functional of the whole system [^9[ . 

At high temperatures the minimum of L(M, 0) is at 
M = 0. At low temperatures the point M = corre- 
sponds to the maximum of this functional. Let us expand 
L(M, 0) with respect to M 

L(M, 0) = L (T) - L 2 (T)M 2 + O (M 4 ), 

L 2 (T) = n 2 (T)-TS 2 . (26) 

Here the coefficients are defined from the the second- 
order expansions of S(M) 



S(M) = S - S 2 M 2 + . . . , 
where S 2 = |iV @, and of fi(m) 

O(m) = n (T) - 2 (T)m 2 + 



(27) 



(28) 



which is to be constructed. Thus the critical tempera- 
ture Tc, below which the functional has the minimum at 
some M/0, is defined from the equation 



L 2 {T C ) =tt 2 {T c )-T c S 2 



0. 



(29) 



The magnetic susceptibility in the PM phase is given 
by the equation 



N 



X(T) 



TS 2 -Q 2 (T) 



(30) 



The grand canonical potential of the electron subsys- 
tem is given by 



fi(m) 



N 
'tt/3 



+ 00 



In 



1 



N(E)dE, (31) 



where N{E) is the DOS given by Eq. ( |ll| ) and \i is the 
chemical potential determined from the equation 



+ OC 



f(E >f i)N(E)dE, 



(32) 



in which f(E, fi) is the Fermi function and n = 1 — x is 
the number of electrons per site. 



V. CALCULATION OF T c AND X (T) 

Due to the constancy of n (or x) the calculation of 
£l 2 (T) requires only the second-order expansion of N(E) 
with respect to m, which is obtained via Eq. ( [ll|) thus 
expanding Eq. (|l8|) 



(33) 



Substituting the related result for A^( 2 ) (E) into Eqs. Q , 
(p9|) and (|3^) , we obtain the following equations 



T C = 0(T C ^(T C )), 



3 1 

iV X( ' ~ T - S(T, n(T)) 



(34) 



(35) 



In these equations n(T) is determined from Eq. ( |32 
where N(E) is replaced by its PM value N^(E), and 



e(T )M ) = |^ Tin 



Im 



1 + e HJ(E-/0 
(5,)e <fl?>« + ¥ ((^(5 £ 3 ) £ -<9 £ 2 ) 



(1- W ( 5 2 )e) (i_ f {g 2 )e y 



dE, (36) 



where 



1 



9t 



E — e — wg. 



(0) 



(37) 



Using Eq. ( |20| ) and integrating by parts, we obtain 

/oc 
f(E,n)8(E)dE, (38) 
-co 

where 

■)/> I fn.-\ I 

(39) 



(/,:/,:, fl.ni ;,/ ' 



3^ \ 1 - f (5?) e J ' 

The function 9(E) being integrated with respect to en- 
ergy gives exactly zero. This leads to Tc(x) =Tc(l — x), 
which reflects the particle-hole symmetry of our model 
irrespective of the disorder strength. 

In the case where the DOS is smooth, the inequality 
W ^> Tc allows us to consider electrons as nearly degen- 
erate. In this case both /i and 0(T, fx) do not depend 
upon T, and Eq. ( p4| ) is just a ready formula for Tq. 
The same is true if there is developed gap. In this case 
fi is near the middle of the gap, so we can substitute 
f(E, n) by one, provided that the integration in Eq. ( pSj ) 
extends only over the filled sub-band. In both these cases 
x(T) obeys the Curie- Weiss law (see Eq. (|35|)). In the 
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first case the ferromagnetic order is mediated mostly by 
mobile holes that is specific for DE. In the second case 
the concentration of the mobile carriers (both holes and 
electrons) is exponentially small. So effective exchange 
between core spins is mostly due to virtual transitions 
of electrons from the lower filled to the upper empty 
sub-band via the gap. This mechanism is an analog of 
super-exchange (SE) acting in the system with electron 
disorder. 

If A ~ T or there is a pseudogap with a strong dip in 
the DOS, the integration in Eqs. @ and @ should 
take into account the tails of f(E,/j,). The exchange in 
such cases is intermediate between DE and SE types. 

For the case of no on-site disorder our Eq. ( |34|) for 
Tc coincides with Eq. (49) of Ref. ^0|, obtained in the 
framework of DMFA and calculated numerically versus 
x. In this case we even managed to calculate the integral 
in Eq. (B8|) analytically, to get for Tc- 



T c = 



WV2 



4tt 



- 4= tan_1 V3(i - v 2 ) 

V3 



(40) 

y = yfin/W, (41) 



while Eq. 



takes the form 

1 I 

2 ~ 7T 



sin y + y 



(42) 



When the disorder is taken into account the integra- 
tion in Eq. (|38|-(|34]) is done numerically. The results for 
Tc{x) at different V/W are presented on Fig. 4. The 
upper (dashed) curve calculated at V = is the same as 
plotted using Eqs. © and @ 
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FIG. 4. The calculated Curie temperature as function of x 
at different V/W. 

One notices that the increase of the impurity poten- 
tial leads to progressive decrease of Tc- This decrease 
becomes more substantial at V/W and x able to pro- 
duce the gap in the density of states. This trend marks 



the modification of the ferromagnetic exchange mecha- 
nism which accompanies the electron band splitting. As 
in the case of DE alone Tc still increases from zero to a 
maximum value upon increasing x from x = to x = 0.5. 



VI. DISCUSSION 

The DE scenario of the FM transition in manganites 
attracted especial attention after its validity was ques- 
tioned in Ref. ^f|. In the following discussion, the po- 
sition of one side can be expressed by reformulation of 
the rhetorical question by Hubbard p2| : "How, in the 
itinerant model, explain a Curie temperature ~ 1000° 
for iron, when calculations always give an exchange field 
^1 — 2 eV?" in the form: How, in the DE model, ex- 
plain a Curie temperature ~ 300° for a manganite (at 
optimal doping), when calculations always give a band- 
width ^1 — 2 eV? The fallacy is based on the implicit 
assumption that the Tc should be of the order of the en- 
ergy difference between the fully ordered FM state and 
the PM state (which is the smallest of the band width 
and the exchange integral) . It does not take into account 
that the fluctuations which restore magnetic symmetry 
in itinerant models are local fluctuations in the orienta- 
tion of magnetic moments, with much lower energy p2[ . 
That is why the numerical coefficient in Eq. ([if)]) turns 
out to be much less than one. 

Our results show that Tc is further decreased by the 
presence of impurities with strong enough potential. It 
may be questioned in this connection why low value of 
Tc in manganites would be on account of fluctuations 
beyond DMFA in the pure DE model rather than due to 
the alloy disorder. The present study reveals an inter- 
esting, though hardly experimentally detectable, feature 
- the alternation of ferromagnetic exchange mechanism 
from DE to SE like as the impurity band splits off the 
conduction band. 
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